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Abstract. Prony mapping provides the global solution of the Prony system 
of equations 

E? =1 A<xJ = m k , k = 0, 1, . . . , 2n - 1. 
This system appears in numerous theoretical and applied problems arising in 
Signal Reconstruction. The simplest example is the problem of reconstruction 
of linear combination of <5-functions of the form g(x) = ^27—1 a i^( x ~ 
with the unknown parameters a;, Xi, i = 1, . . . , n, from the "moment mea- 
surements" m,fe = J x k g(x)dx. 

Global solution of the Prony system, i.e. inversion of the Prony mapping, 
encounters several types of singularities. One of the most important ones is a 
collision of some of the points . The investigation of this type of singularities 
has been started in 1211 where the role of finite differences was demonstrated. 

In the present paper we study this and other types of singularities of the 
Prony mapping, and describe its global geometry. We show, in particular, close 
connections of the Prony mapping with the "Vieta mapping" expressing the 
coefficients of a polynomial through its roots, and with hyperbolic polynomials 
and "Vandermonde mapping" studied by V. Arnold. 



1. Introduction 

Prony system appears as we try to solve a very simple "algebraic signal reconstruc- 
tion" problem of the following form: assume that the signal F(x) is known to be a 
linear combination of shifted <5-functions: 

d 

F(x)=J2a j 6(x-x j ). (1.1) 

i=i 

We shall use as measurements the polynomial moments: 

m k =m k {F) = f x k F(x)dx. (1.2) 



After substituting F into the integral defining m/. we get 

d d 



mfc(F) = / x k ajS(x — Xj) d x = ajXj. 
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Considering Oj and unknowns, we obtain equations 

d 

m k (F) =^a j x k j , A; = 0,1,.... (1.3) 
j'=i 

This infinite set of equations (or its part, for k = 0, 1, . . . , 2d — 1), is called Prony 
system. It can be traced at least to R. de Prony (1795, [T5]) and it is used in a 
wide variety of theoretical and applied fields. See [2] for an extensive bibligoraphy 
on the Prony method. 

In writing Prony system (jl.3[) we have assumed that all the nodes x\,...,Xd are 
pairwise different. However, as a right-hand side fi = (mo, . . . , 7Ti2d-i) of ([1.3)1 is 
provided by the actual measurements of the signal F, we cannot guarantee a priori, 
that this condition is satisfied for the solution. Moreover, we shall see below that 
multiple nodes may naturally appear in the solution process. In order to incorporate 
possible collisions of the nodes, we consider "confluent Prony systems". 

Assume that the signal F(x) is a linear combination of shifted <5-functions and their 
derivatives: 

FW^^a/)^-^). (1.4) 

Definition 1.1. For F (x) as above, the vector D (F) = (d\, . . . , d s ) is the multi- 
plicity vector of F, s — s (F) is the size of its support, T (F) = f (xi, . . . , x s ), and 
rank (F) = f X)j=i dj is its rank. For avoiding ambiguity in these definitions, it is 
always understood that a^^ .-i 7^ for all j = 1, . . . , s (i.e. dj is the maximal index 
for which aj^—i 7^ 0). 

For the moments ra^ = rrik(F) = J x k F{x) dx we now get 

3=1 1=0 v '' 
Considering Xi and ajj as unknowns, we obtain a system of equations 

E E 7fc37fT a ^4"" = m *> k = 0,l,...,2d-l, (1.5) 
3=1 £=0 ^ 

which is called a confluent Prony system of order d with the multiplicity vector 
D = (di, . . . , d s ). The original Prony system (| 1 . 3[) is a special case of the confluent 
one, with D being the vector (1, . . . , 1) of length d. 

The system (|1.5|) arises also in the problem of reconstructing a planar polygon P 
(or even an arbitrary semi-analytic quadrature domain) from its moments 

m-kixp) = // -2 fe Xpdxdy, z = x + ly, 
JJm? 

where xp is the characteristic function of the domain P C M. 2 . This problem is 
important in many areas of science and engineering [11] . The above yields the 
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confluent Prony system 

s dj-l 

m k = 53 c ^ k ( k -!)•••(*-* + 1)2*"*, Cij eC, Zj eC\ {0} . 

j=l i=0 

Definition 1.2. For a given multiplicity vector D = (d\, . . . , d s ), its order is 

As we shall see below, if we start with the measurements fi(F) = fj, = (mo, . . . , m 2( i_i), 
then a natural setting of the problem of solving the Prony system is the following: 

Problem 1.3 (Prony problem of order d). Given the measurements 

(i = (mo, ■ • ■ , m 2 d-i) G C 2d 
m fte right hand side of (|1.5p . /ind i/ie multiplicity vector D = (dx, . . . , d s ) of order 
r = ^2j=i dj < d, and find the unknowns Xj and which solve the corresponding 
confluent Prony system (|1.5p with the multiplicity vector D (hence, with solution 
of rank r ). 

It is extremely important in practice to have a stable method of inversion. Many 
research efforts are devoted to this task (see e.g. [31 [71 EH1 HH 120] and references 
therein). A basic question here is the following. 

Problem 1.4 (Noisy Prony problem). Given the noisy measurements 

fi = (m , ■ • ■ ,m 2 d-i) G C 2d 

and an estimate of the error \ fhk — rnk\ < £fe, solve Problem \1.3\ so as to minimize 
the reconstruction error. 

In this paper we study the global setting of the Prony problem, stressing its alge- 
braic structure. In Section [5] the space where the solution is to be found (Prony 
space) is described. It turns out to be a vector bundle over the space of the nodes 
x\, . . . , Xd- We define also three mappings: "Prony", "Taylor", and "Stieltjes" ones, 
which capture the essential features of the Prony problem and of its solution pro- 
cess. 

In Section[3]we investigate solvability conditions for the Prony problem. The answer 
leads naturally to a stratification of the space of the right-hand sides, according to 
the rank of the associated Hankel-type matrix and its minors. The behavior of the 
solutions near various strata turns out to be highly nontrivial, and we present some 
initial results in the description of the corresponding singularities. 
In SectionSl we study the multiplicity-restricted Prony problem, fixing the collision 
pattern of the solution, and derive simple bounds for the stability of the solution 
via factorization of the Jacobian determinant of the corresponding Prony map. 
In Section [5] we consider the rank-restricted Prony problem, effectively reducing 
the dimension to 2r instead of 2d, where r is precisely the rank of the associated 
Hankel-type matrix. In this formulation, the Prony problem is solvable in a small 
neighborhood of the exact measurement vector. 

In Section[6]we study one of the most important singularities in the Prony problem: 
collision of some of the points Xi . The investigation of this type of singularities has 
been started in [21] where the role of finite differences was demonstrated. In the 
present paper we introduce global bases of finite differences, study their properties, 
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and prove that using such bases we can resolve in a robust way at least the linear 
part of the Prony problem at and near colliding configurations of the nodes. 

In Section [7] we discuss close connections of the Prony problem with hyperbolic 
polynomials and "Vandermonde mapping" studied by V.I. Arnold in pQ and by 
V.P.Kostov in PJ5J HH [H], and with "Vieta mapping" expressing the coefficients 
of a polynomial through its roots. We believe that questions arising in theoretical 
study of Prony problem and in its practical applications justify further investigation 
of these connections, as well as further applications of Singularity Theory. 

Finally, in Appendix [A] we describe a solution method for the Prony system based 
on Pade approximation. 

2. Prony, Stieltjes and Taylor Mappings 

In this section we define "Prony", "Taylor", and "Stieltjes" mappings, which capture 
some essential features of the Prony problem and of its solution process. The main 
idea behind the spaces and mappings introduced in this section is the following: 
associate to the signal F{x) — J^i=i a%S(x — xi) the rational function R(z) — 
J2i=i g-x • f ac t' ^ S ^ ne Stieltjes integral transform of F). The functions 
R obtained in this way can be written as R(z) — with degP < degQ — 1, 
and they satisfy R(oo) = 0. Write R as R(z) = J2i=i i- "'/z • Developing the 
summands into geometric progressions we conclude that R(z) — Y^kLo m k(^;) k+1 , 
with mfc = ^2f =1 ciiXi, so the moment measurements in the right hand side 
of the Prony system (jl.3[) are exactly the Taylor coefficients of R(z). We shall 
see below that this correspondence reduces solution of the Prony system to an 
appropriate Pade approximation problem. 

Definition 2.1. For each w = (xi, . . . , xj) <E C d , let s = s (w) be the number 
of distinct coordinates Tj, j — l,...,s, and denote T (w) — (ti,...,t s ). The 
multiplicity vector is D = D (w) = (di, ■ ■ ■ , d s ), where dj is the number of times 
the value Tj appears in {xi, . . . , Xd} ■ The order of the values in T (w) is defined by 
their order of appearance in w. 

Example 2.2. For w = (3, 1, 2, 1, 0, 3, 2) we have s (w) = 4, T (w) = (3, 1, 2, 0) and 
D{w) = (2,2,2,1). 

Remark 2.3. Note the slight abuse of notations between Definition 11.11 and Defini- 
tion [2Tj Note also that the order of D (w) equals to d for all w € C d . 

Definition 2.4. For each w € C d , let s = s (w) , T (w) = (n, . . . , r s ) and D (w) — 
(di, . . . , d s ) be as in Definition 12. II 

(1) V w is the vector space of dimension d containing the linear combinations 

s dj — 1 

of (5-functions and their derivatives at the points of T(w). The "standard 
basis" of V w is given by the distributions 

8 i 4 = 6F>{x-T j ), j = l,...,s(w); i = 0,...,dj-l. (2.2) 
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(2) W w is the vector space of dimension d of all the rational functions with 
poles T (w) and multiplicities D (w) , vanishing at oo : 

R ^ = tM' Q(*) = i\{*-Ti) di > degP(z)<degQ^d. 
Q ( z ) fj[ 

The "standard basis" of W w is given by the elementary fractions 
Rj,e = - — —-j, j = l,...,s; £ = l,...,dj. 

(Z - Tj) 

Now we are ready to formally define the Prony space Vd and the Stieltjes space Sd- 

Definition 2.5. The Prony space Vd is the vector bundle over C d , consisting of 
all the pairs 

(w,g) : w e C d , g G V w . 

The topology on Vd is induced by the natural embedding Vd C C d x V, where V is 
the space of distributions on C with its standard topology. 



Definition 2.6. The Stieltjes space Sd is the vector bundle over C d , consisting of 
all the pairs 

(to, 7): weC d ,-ieW w . 

The topology on Sd is induced by the natural embedding Sd C C d x 1Z, where 1Z is 
the space of complex rational functions with its standard topology. 



Definition 2.7. The Stieltjes mapping SM : Vd — > Sd is defined by the Stieltjes 
integral transform: for (w,g) € Vd 

g(x)dx 



f 

SM((w,g)) = (u>,7), l(z)= / 



Sometimes we abuse notation and write for short SM (g) = 7, with the under- 
standing that SM. is also a map SM :V W ->W W for each w £ C d . 

The following fact is immediate consequence of the above definitions. 

Proposition 2.8. SM is a linear isomorphism of the bundles Vd and Sd (for each 
w G C d , SM is a linear isomorphism of the vector spaces V w and W w ). In the 
standard bases ofV w and W w , the map SM is diagonal, satisfying 

sm (5 jtt ) = {-ifem^ (z). 

Furthermore, for any (w, g) G Vd 

SM{g)= ^4 ' degP<dcgQ = rank(.g) < d. (2.3) 

irreducible 

Definition 2.9. The Taylor space Td is the space of complex Taylor polynomials 
at infinity of degree 2d — 1 of the form Yl'k^o m k{\) k+1 - We shall identify Td with 
the complex space C 2d with the coordinates mo, . ■ . , mid-i- 
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Definition 2.10. The Taylor mapping TM : Sd — > Td is defined by the truncated 
Taylor development at infinity: 

2d-i / \\ k+1 00 /l\ fe+1 

TM ((w, 7)) = ^ q/g ( - J where 7 (z) = ^ a k y -J 



,z / * — ' Vz 

fc=0 fc=0 



We identify T.M ((w, 7)) as above with (ao, ■ • ■ , a^d-i) € C 2d . Sometimes we write 
for short TM (7) = (a , . . . , a 2 d-i)- 

Finally, we define the Prony mapping VM which encodes the Prony problem. 

Definition 2.11. The Prony mapping VM : Vd C 2d for (w,g) G is defined 
as follows: 

VM ((w.g)) = (m , . . . ,m 2d -i) G C 2d , m k = m k (g) = Jx k g(x)dx. 

By the above definitions, we have 

VM = TMoSM. (2.4) 

Solving the Prony problem for a given right-hand side (too, ■ • ■ ; i r n-2d-i) is therefore 
equivalent to inverting the Prony mapping VM. As we shall elaborate in the sub- 
sequent section, the identity (|2.4|) allows us to split this problem into two parts: 
inversion of TM, which is, essentially, the Pade approximation problem, and in- 
version of SM, which is, essentially, the decomposition of a given rational function 
into the sum of elementary fractions. 



3. Solvability of the Prony problem 

3.1. General condition for solvability. In this section we provde a necessary 
and sufficient condition for the Prony problem to have a solution (which is unique, 
as it turns out by Proposition I3.2p . As mentioned in the end of the previous 
section, our method is based on inverting (I2.4[) and thus relies on the solution of 
the corresponding (diagonal) Pade approximation problem [1]. 

Problem 3.1 (Diagonal Pade approximation problem). Given /i = (mo, . . . , m^d-i) G 
C 2d , find a rational function Rd{z) = G Sd with dcgP < degQ ^5 d, such that 
the first 2d Taylor coefficients at infinity of Rd(z) are {mfc}?^ 1 . 

Proposition 3.2. A solution to Problem \S.l[ if exists, is unique. 

Proof. Writing R(z) = 04, i?i (z) = with deg P < deg Q ^ d and deg Pi < 

degQi < d, we get 

R-R x = PQl - PlQ 



QQi ' 

and this function, if nonzero, can have a zero of order at most 2d — I at infinity. □ 



Let us summarize the above discussion with the following statement. 



GEOMETRY AND SINGULARITIES OF THE PRONY MAPPING 



7 



Proposition 3.3. The tuple 

js, D = (d u ...,d s ), r = Y^dj<d, X = {xj} a j=1 , A {».,., | , , s . , u d : | 

is a ( unique, up to a permutation of the nodes {xj } ) solution to Problem \1.3\ with 
right-hand side 

H = (mo, • ■ • ,m 2 d_i) e C 2d 
if and only if the rational function 

s dj 2d-l 
j=l fcl l Z x j) k=0 

is a (unique) solution to Problem \3.1\ with input fi. In that case, 

/oo (x)dx s dj ~ 1 
where g (x) = V] Y] (a; - aij) , 

i.e. Rd,x.A (z) is the Stieltjes transform of g{x). 

Proof. This follows from the definitions of Section [51 (|2.4[) , Proposition 13.21 and 
the fact that the problem of representing a given rational function as a sum of 
elementary fractions of the specified form (i.e. inverting SAA) is always uniquely 
solvable up to a permutation of the poles. □ 

The next result provides necessary and sufficient conditions for the solvability of 
Problem 13.11 It summarizes some well-known facts in the theory of Pade approx- 
imation, related to "normal indices" (see, for instance, [1]). However, these facts 
are not usually formulated in the literature on Pade approximation in the form we 
need in relation to the Prony problem. Consequently, we give a detailed proof of 
this result in Appendix [XJ This proof contains, in particular, some facts which are 
important for understanding the solvability issues of the Prony problem. 

Definition 3.4. Given a vector /i = (mo, . . . , to^-i), let Md denote the d X (d + 1) 
Hankel matrix 



Mr, 



to toi m 2 ... m d 
TOi m 2 to 3 ... rrid+i 

TO d _i m d m d+ i . . . m 2 d-i 



(3.1) 



For each e ^ d, denote by M e the ex (e + 1) submatrix of Md formed by the first 
e rows and e + 1 columns, and let M e denote the corresponding square matrix. 

Theorem 3.5. Let fi = (mo, . . . , iJi2d-i) be given, and let r ^ d be the rank of the 
Hankel matrix Md as in (|3.ip . Then Problem \3.1\ is solvable for the input fi if and 
only if the upper left minor \M r \ of Md is non-zero. 



As an immediate consequence of Theorem 13.51 and Proposition 13.31 we obtain the 
following result. 
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Theorem 3.6. Let fi = (too, • . . , m,2d-i) be given, and let r ^ d be the rank of 
the Hankel matrix Md as in (|3.1jl . Then Problem \1.3\ with input fj, is solvable if 
and only if the upper left minor \M r \ of Md is non-zero. The solution, if exists, 
is unique, up to a permutation of the nodes {xj}. The multiplicity vector D = 
(di, . . . , d s ), of order X)?=i dj = T, of the resulting confluent Prony system of rank 
r is the multiplicity vector of the poles of the rational function Rd,x,A (z), solving 
the corresponding Pade problem. 

As a corollary we get a complete description of the right-hand side data fi 6 C 2d 
for which the Prony problem is solvable (unsolvable). Define for r = 1, . . . , d sets 
E r C C 2d (respectively, Y' r C C 2d ) consisting of /i £ <C 2d for which the rank 
of Md — r and \M r \ ^ (respectively, \M r \ = 0). The set E r is a difference 
E r = E* \ E 2 of two algebraic sets: E* is dehned by vanishing of all the sx s minors 
of Md, r < s < d, while E 2 is defined by vanishing of \M r \. In turn, Y' r = E r x \ E r 2 , 
with E,. 1 = n E 2 and E r 2 defined by vanishing of all the r x r minors of Md- The 
union E r U Y' r consists of all /i for which the rank of Md = r, which is \ E r 2 . 

Corollary 3.7. The set E (respectively, Y! ) of n € C 2d for which the Prony problem 
is solvable (respectively, unsolvable) is the union E = uf =1 Y r (respectively, E' = 
U^ =1 S; ). In particular, E' C {/i G C 2d , det M d = 0}. 

So for a generic right hand side /i we have |M<j| ^ 0, and the Prony problem 
is solvable. On the algebraic hypersurface of fi for which |M<j| = 0, the Prony 
problem is solvable if Md-i ^ 0, etc. 

Let us now consider some examples. 

Example 3.8. Let us fix d = 1,2,.... Consider /i = (mo, . . . , TO2d-i) € C 2d , 
the right hand sides of the Prony problem, to be of the form fj, = m = (6m) 
Mi i ) . , 0, 0), with all the nik = besides mi — 1, t= 0, 2d— 1, 

position £+1 

and let Aij be the corresponding matrix. 

Proposition 3.9. The rank of Mj is egwaZ to £ + 1 for I < d— 1, and if is eouaZ 
to 2d — £ for £ > d. The corresponding Prony problem is solvable for £ < d — 1, and 
it is unsolvable for £ > d. 

Proof. For d — 5 and ^ = 2, 4, 5, 9, the corresponding matrices M d are as follows. 



Ml 



Ml 









1 













"0 











1 








1 

























1 








1 

















. Mi = 








1 


































1 





























0_ 




1 
































f 




"0 





























1 


































1 


































1 


































1 































1 



(solvable) 



(unsolvable) 
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In general, the matrices M d have the same pattern as in the special cases above, so 
their rank is I + 1 for £ *C d — 1, and 2d — £ for £ ^ d, as stated above. Application 
of Theorem 13.61 completes the proof. □ 

In fact, fi£ is a moment sequence of 

F(x) = ±5W (x), 

and this signal belongs to Vd if and only if I ^ d — 1. In notations of Corollary 13. 71 
we have 

Hi e T, e+ i, l < d - 1, 

It is easy to provide various modifications of the above example. In particular, for 
fj, = fn = (0,...,0,1,1,...,1), the result of Proposition 13.91 remains verbally true. 

Example 3.10. Another example is provided by fj,t lt t 3 , with all the mu — besides 
m ti = m t2 = 1, < t\ < d < 4 < 2d - 1. For l\ < £ 2 - d + 1 the rank of 
the correspondent matrix Md is r = 2d + £\ — £2 + 1 while \M r \ = 0, so the Prony 
problem for such fJH lt e. 2 is unsolvable. For d — 5 and £\ = 2, £2 = 8 the matrix is 
as follows: 

0" 

. 

1 

1 0_ 

3.2. Near-singular inversion. The behavior of the inversion of the Prony map- 
ping near the unsolvability stratum £' and near the strata where the rank of Md 
drops, turns out to be pretty complicated. In particular, in the first case at least 
one of the nodes tends to infinity. In the second case, depending on the way the 
right-hand side /1 approaches the lower rank strata, the nodes may remain bounded, 
or some of them may tend to infinity. In this section we provide one initial result 
in this direction, as well as some examples. A comprehensive description of the in- 
version of the Prony mapping near E' and near the lower rank strata is important 
both in theoretical study and in applications of Prony-like systems, and we plan to 
provide further results in this direction separately. 

Theorem 3.11. As the right-hand side /1 6 C 2d \ S' approaches a finite point 
/jq 6 S'| at least one of the nodes X\ , . . . , Xd in the solution tends to infinity. 

Proof. By assumptions, the components mo, . . . , rri2d~i of the right-hand side fi = 
(mo, . . . ,m2d-i) € C 2d remain bounded as fi — > fio- By Theorem 16.141 the finite 
differences coordinates of the solution VA4~ 1 (n) remain bounded as well. Now, if 
all the nodes are also bounded, by compactness we conclude that VMT — > oj g 
Vd- By continuity in the distribution space fLemma 16. 6[) we have VM.(iS) = /io- 
Hence the Prony problem with the right-hand side /io has a solution lu 6 Vd, in 
contradiction with the assumption that /jq 6 S'. □ 









1 








1 








1 
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Example 3.12. Let us consider an example: d — 2 and fio = (0, 0, 1, 0). Here the 
rank t of M 2 is 2, and \M 2 \ = 0, so by Theorem l3.6l we have /j,q £ T! 2 C Consider 
now a perturbation fi(e) = (0, e, 1, 0) of fiQ. For e^Owe have fi(e) 6 S2 C S, and 
the Prony system is solvable for /i c . Let us write an explicit solution: the coefficients 
Co, ci of the polynomial Q{z) = cq + c\z + z 2 we find from the system (|AI**|k 



"0 e" 




co 






e 1. 




_Cl. 




_ . 



whose solution is c\ = —-, cq = Hence the denominator Q(z) of R(z) is 
Q(z) = — iz + z 2 , and its roots are x\ — 1+ 2 , X2 = ^ ■ The coefficients 
60, ^1 of the numerator P(z) = 6q + &i z we hud from (|A1*[) : 



i.e. bi 



-0, b Q 
R(z)-- 



e. Thus the solution of the associated Pade problem is 











-_v 

e 









( I. 




_ 1 . 




M 



P{z) 



1 



1 



Q(z) (z - xi)(z - x 2 ) iy/3{z-xi) i^{z-x 2 )' 
Finally, the (unique up to a permutation) solution of the Prony problem for /i e is 



a 1 



e 

V3' 



a 2 = - - 



Xl = 



1 + zV3 



X2 = 



As e tends to zero, the nodes x\, x 2 tend to infinity while the coefficients a\, a 2 tend 
to zero. 

As it was shown above, for a given /i € S (say, with pairwise different nodes) the 
rank of the matrix Md is equal to the number of the nodes in the solution for which 
the corresponding (^-function enters with a non-zero coefficients. So fi approaches 
a certain /j,q belonging to a stratum of a lower rank of Md if and only if some of 
the coefficients aj in the solution tend to zero. We do not analyze all the possible 
scenarios of such a degeneration, noticing just that if /zo G S', i.e., the Prony 
problem is unsolvable for /j,q, then Theorem 13. 1 II remains true, with essentially the 
same proof. So at least one of the nodes, say, Xj, escapes to infinity. Moreover, 
one can show that a,jX 2d ~ Y cannot tend to zero - otherwise the remaining linear 
combination of 8- functions would provide a solution for /iq. 

If /io £ S, i.e., the Prony problem is solvable for /io, all the nodes may remain 
bounded, or some Xj may escape to infinity, but in such a way that ajX 2 ^ 1 tends 



to zero. 



4. Multiplicity-restricted Prony problem 

Consider Problem 11.41 at some point fiQ € S. By definition, /io <= S ro for some 
ro < d. Let ^0 — VAi ((wo,5o)) f° r some (wo,ffo) S Td- Assume for a moment that 
the multiplicity vector Dq = D (go) = (d±, . . . d So ), 52^==i ^1 ~ r °' ^ a3S a non_ trivial 
collision pattern, i.e. dj > 1 for at least one j = 1, . . . , 8q. It means, in turn, that 
the function Rd ,x.a {z) has a pole of multiplicity dj. Evidently, there exists an 
arbitrarily small perturbation jl of fio for which this multiple pole becomes a cluster 
of single poles, thereby changing the multiplicity vector to some D' ^ Do. While 
we address this problem in Section [6] via the bases of divided differences, in this 
section we consider a "multiplicity-restricted" Prony problem. 
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Definition 4.1. Let x = (xi, ...,i s )eC s and D — (di, . . . , d s ) with d = Ylj=i dj 
be given. The d x d confluent Vandermonde matrix is 



V = V(x,D) = V(x 1 ,d 1 ,...,x s ,d s ) 



Vl,0 V2,0 ■ • • v Sj o 
Vl,l v 2 ,i ... v M 

Vl,d-1 V 2 ,d-1 • • • V s d-1 



(4.1) 



where the symbol vj k denotes the following 1 x dj row vector 

dcf f u , fc_1 , /, ,n /, i \ k-dj+1 



v .i.k 



[ x$, kx k ~\ ... ,k(k-i)---(k-d j )x k r dj+1 ] . 



Proposition 4.2. The matrix V defines the linear part of the confluent Prony 
system (jl.5[) in the standard basis for V w , namely, 



(4.2) 



Definition 4.3. Let VM. (wo,go) = Ho €E S ro with D (go) = Do and s (go) = So- 
Let Vd denote the following subbundle of Vd of dimension so + ro' 

V Do ={(w,g) eV d : D(g)=D Q }. 

The multiplicity-restricted Prony mapping VA4* Do : Vd ~^ C 50 "*" 1 * is the composi- 
tion 

VM* Do =woVM \v Do , 
ya+r a j s (- ne projection map on the first sq + r o coordinates. 





Ol,0 




to 








TOl 


V (xi,di, . . .,x a ,d a ) 










_a s ,d s -l_ 




TOd-1. 



where 7r : C 



2</ 



Inverting this T j M.* Dq represents the solution of the confluent Prony system (ll.5[) 
with fixed structure Do from the first k — 0, 1, . . . , so + r o ~ 1 measurements. 

Theorem 4.4 ([7]). Lei = VM* D() ((wo,go)) € C So+rci wif/i i/je unperturbed solu- 
tion go — ^2j=i Sfco 1 a j,#fi^ ( x ~ T j)- I n a sma M neighborhood of (wo,go) € *Pd , 
the map VM* Dg is invertible. Consequently, for small enough e, the multiplicity- 
restricted Prony problem with input data fx* € t£ r o+so satisfying \\fl* — fj,o\\ < e has 
a unique solution. The error in this solution satisfies 



|Aa,- 



I At, 



< 



< 



2 
J\ 

2 



so+ro 



so+ra 



SO + r a 



l a J,l-l| 



\ a 3,d 3 



where 5 



def 



Tj| (for consistency we take Oj.-i — in the above formula) . 



Proof outline. The Jacobian of V M* Dq can be easily computed, and it turns out to 
be equal to the product 



Jvm- =V(Ti,di + l, 



,r So ,d So + l)diag{J5,-} 
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where V is the confluent Vandermonde matrix (|4.1| on the nodes (n , . . . , r So ) , with 
multiplicity vector 

£>o = (d 1 + l,...,d, +l), 
while -E is the (dj + 1) x (dj + 1) block 

"1 ••• 
1 ••• aj-o 







Since /iq G S r , the highest order coefficients a^dj-x are nonzero. Furthermore, since 
all the Tj are distinct, the matrix V is nonsingular. Local invertability follows. To 
estimate the norm of the inverse, use bounds from [5J. □ 

Remark 4.5. Note that as two nodes collide (S — > 0), the inversion of the multiplicity- 
restricted Prony mapping VM* Do becomes ill-conditioned proportionally to <5 _ ( s °+ r °) . 

Let us stress that we are not aware of any general method of inverting 'PA4* Do , 
i.e. solving the multiplicity-restricted confluent Prony problem with the smallest 
possible number of measurements. As we demonstrate in [5], such a method exists 
for a very special case of a single point, i.e. s = 1. 



5. Rank-restricted Prony problem 

Recall that the Prony problem consists in inverting the Prony mapping VA4 : Vd — > 
Td- So, given y, = (mo, . . . ,W2d-i) £ Td we are looking for {w,g) £ Vd such that 
m k{g) — j x k g(x)dx — nik, with k = 0, 1, . . . , 2d — 1. If /i € £ r with r < d, then 
in fact any neighborhood of /i will contain points from the non-solvability set £'. 
Indeed, consider the following example. 



Example 5.1. Slightly modifying the construction of Example l3.101 consider fig lt i 2 ^ 

1 and mi 2 = e, such that £2 > l\ + d — 1. 
, the corresponding matrix is 



C 2d with all the mk = besides me 1 
For example, if d = 5 and £\ = 2, £2 



M 5 (2 ' 8 '' 



For e = the Prony problem is solvable, while for any small perturbation e ^ 
it becomes unsolvable. However, if we restrict the whole problem just to d = 3, it 
remains solvable for any small perturbation of the input. 

We therefore propose to consider the rank-restricted Prony problem analogous to 
the construction of Section |4l but instead of fixing the multiplicity D (g) we now 
fix the rank r (recall Definition II. 1| . 

Definition 5.2. Denote by V r the following vector bundle: 

V r = {(w,g): toeC, geV w }, 
where V w is defined exactly as in Definition 12.41 replacing d with r. 



Likewise, we define the Stieltjes bundle of order r as follows. 
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Definition 5.3. Denote by S r the following vector bundle: 

where W w is defined exactly as in Definition 12.41 replacing d with r. 



The Stieltjes mapping acts naturally as a map SM : V r — > S r with exactly the 
same definition as Definition 12.71 

The restricted Taylor mapping TM r ■ S r — > C 2r is, as before, given by the trun- 
cated development at infinity to the first 2r Taylor coefficients. 

Definition 5.4. Let tt : C 2d — > C 2r denote the projection operator onto the first 
2r coordinates. Denote S* = f 7r (S r ). The rank-restricted Prony mapping V : M* : 
V r —> £* is given by by 

VM* {(w,g)) = (mo, ■ ■ ■ ,m 2r -i) , m k = m k (g) = Jx k g(x)dx. 

Remark 5.5. V r can be embedded in Vd, for example by the map S r : V r — > Vd 

E r : (w,g) e V r i — >■ {w',g') eV d - w' = I xi, . . . ,x r , 0, . . .0 J , g' = p. 

\ x(d-r)/ 

With this definition, PA^* can be represented also as the composition 

VM*. =7roPMoS r . 
Proposition 5.6. TTie rank-restricted Prony mapping satisfies 

VM* = TM r oSM. 

Inverting VM.* represents the solution of the rank-restricted Prony problem. Unlike 
in the multiplicity-restricted setting of Section 31 here we allow two or more nodes 
to collide (thereby changing the multiplicty vector D (g) of the solution). 

The basic fact which makes this formulation useful is the following result. 

Theorem 5.7. Let /ij € S*. Then in a small neighborhood of ^ £ C 2r , the Taylor 
mapping T M r is continuously invertible. 

Proof. This is a direct consequence of the solution method to the Pade approxi- 
mation problem described in Appendix |A] Indeed, if the rank of M r is full, then it 
remains so in a small neighborhood of the entire space C 2r . Therefore, the system 
(|Al*-*r[) remains continuously invertible, producing the coefficients of the denomina- 
tor Q (z). Consequently, the right-hand side of (IA1*I) depends continuously on the 
moment vector n* = (mo, . . . , m^r-i) € C 2r . Again, since the rank always remains 
full, the polynomials P (z) and Q (z) cannot have common roots, and thereby the 
solution R = ^ = TMr 1 (fi*) depends continuously on /i* (in the topology of the 
space of rational functions) . □ 



In the next section, we consider the remaining problem: how to invert SM in this 
setting. 
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6. Collision singularities and bases of finite differences 

6.1. Introduction. Collision singularities occur in Prony systems as some of the 
nodes Xi in the signal F(x) = Yli=i a i$( x ~ x i) approach one another. This happens 
for /i near the discriminant stratum A C C 2d consisting of those (mo, ■ ■ ■ , W2d-i) 
for which some of the coordinates {xj} in the solution collide, i.e. the function 
Rd.x,a (z) has multiple poles (or, nontrivial multiplicity vector D). As we shall 
see below, typically, as (i approaches /iq G A, i.e. some of the nodes Xi collide, the 
corresponding coefficients tend to infinity. Notice, that all the moments rrik = 
mk{F) remain bounded. This behavior creates serious difficulties in solving "near- 
colliding" Prony systems, both in theoretical and practical settings. Especially 
demanding problems arise in the presence of noise. The problem of improvement of 
resolution in reconstruction of colliding nodes from noisy measurements appears in 
a wide range of applications. It is usually called a "super-resolution problem" and 
a lot of recent publications are devoted to its investigation in various mathematical 
and applied settings. See [5] and references therein for a very partial sample. 
Here we continue our study of collision singularities in Prony systems, started in 
|21j . Our approach uses bases of finite differences in the Prony space V r in order 
to "resolve" the linear part of collision singularities. In these bases the coefficients 
do not blow up any more, even as some of the nodes collide. 

Example 6.1. Let r = 2, and consider the signal F — a\5 (x — x\) + 0,26 (x — X2) 
with 

x\ = t, X2 = t + e, 
a\ = -e _1 , a 2 = e _1 . 
The corresponding Prony system is 



{a\x\ + a 2 X2 —) mk 



kth ~ 1+ J2 (j)* fe ~ V_1 ' * = 0,1,2,3. 



dcf 



As e — > 0, the Prony system as above becomes ill-conditioned and the coefficients 
{a.,} blow up, while the measurements remain bounded. Note that 

"0 1 2t + p 2 (t, e) " 

.1 2t + p 2 (t, e) 3t 2 + p 3 (t, e)J ' 



Mn = 



therefore rankM 2 = 2 and \M 2 \ = 1^0, i.e. the Prony problem with input 
(too, ■ • ■ j m 3) remains solvable for all e. However, the standard basis {5 (x — x\) , 5 (x 
degenerates, and in the limit it is no more a basis. If we represent the solution 



in the basis 



F e (x) = --S (x-t) + -6(x-t-e) 
e e 



Ai (£1,2:2) = S(x-xi), 

A 2 (xi,x2) = 5(x — xi)-\ 8 (x — X2) , 

Xi - X 2 X 2 - Xi 



then we have 



F e (x) = l- A 2 (t, t + e), 
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i.e. the coefficients in this new basis are just {61 =0, 61 = 1}. As e — > 0, in fact 
we have 

A 2 (t,t + e)^5' (x-t), 
where the convergence is in the topology of the bundle V r . 

Our goal in this section is to generalize the construction of Example 16.11 and |21| to 
handle the general case of colliding configurations. 

6.2. Divided finite differences. For modern treatment of divided differences, see 
e.g. [HI LLH LIS]- We follow [S] and adopt what has become by now the standard 
definition. 

Definition 6.2. Let an arbitrary sequence of points w — (xi,X2, ■ • ■ , ) be given 
(repetitions are allowed). The (n-l)-st divided difference A™ -1 (w) : II — > C is 
the linear functional on the space II of polynomials, associating to each p E II its 
(uniquely defined) n-th coefficient in the Newton form 

00 i 

P( x ) = ^2 {A^ 1 (xi, . . . ,Xj)p} ■ qj-i, w (x) , q iiW (x) = f (x - xj) . (6.1) 

j=i j=i 

It turns out that this definition can be extended to all sufficiently smooth functions 
for which the interpolation problem is well-defined. 

Definition 6.3 (|9|). For any smooth enough function /, defined at least on 
X\, . . . , x n , the divided finite difference A™ -1 (x\, . . . , x n ) f is the n-th coefficient 
in the Newton form (|6.1|) of the Hermite interpolation polynomial P n , which agrees 
with / and its derivatives of appropriate order on x±, . . . , x n : 

/W ( Xj ) = (x 3 ): 1 < j < n, < I < d j = #{1: x t = Xj } . (6.2) 



Therefore, each divided difference can be naturally associated with an element of 
the Prony space (see Item [S] in Proposition 16.41 and Definition 16.51 below for an 
accurate statement). 

Let us now summarize relevant properties of the functional A which we shall use 
later on. 



Proposition 6.4. For w — (xi, . . . , x n ) £ C n , let s (w) , T (w) and D (w) be 

defir, 

HO 



defined according to Definition \2.1\ Let q r ^ w (z) — Y[j=i ( z — r i) i be defined as in 



(1) The functional A n_1 [x\, . . . , x n ) is a symmetric function of its arguments, 
i.e. it depends only on the set {x%, . . . ,x n } but not on its ordering. 

(2) A™^ 1 (xi, . . . , x n ) is a continuous function of the vector (x\, . . . , x n ). In 
particular, for any test function f 

lim A"" 1 ( Xl ,...,x n )f= A"" 1 (t!,...,t n )/. 

{Xl,—,Xn)—>{tl, — ,t n ) 

(3) A may be computed by the recursive rule 

X\ ^= x ri 




A n ^(x 1 ,...,x n )f={ ( n _ ^ , ^ (6.3) 

3^ A (e,x 2 ,...,x„_i)/} | 5=x „, X!=X n , 

where A (xi) f — f (xi) . 
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(4) Let f z (x) = (z — x) 1 . Then for all z £ {x\ } . . . ,x n } 

A n_1 (xi, x n ) f z = - — —T\- (6-4) 

(5) By (|6.2p . A"^ 1 (xi, . . . ,x n ) is a linear combination of the functionals 

(x - Tj) , 1 < j < s, < I < dj. 
In fact, using (16. 4[) we obtain the Chakalov's expansion (see [9] ) 

s dj-1 

i=i e=o 

where the coefficients {dj/} are defined by the partial fraction decomposi- 
tion^ 

_j_ = la 
Qn, W (z) j^j^ (z- Tj ) e+1 

(6) By (JHIIJ and 

A"' 1 f t^t j = ^-l—^-D (x - t) . (6.7) 

\ xn / 

(7) Popoviciu's refinement lemma [9l Proposition 23]: /or every index subse- 
quence 

1 s$ a (1) < o- (2) < • • • < a (k) < n, 
there exist coefficients a (j) such that 

cr(fc)-fc 

(Xa(l),---,X a (k)) = a ^ ( x 3 + l, x j+2,---,X J+ k)- (6.8) 

J'=°-(1)-1 

Based on the above, we may now identify A with elements of the bundle V r . 

Definition 6.5. Let w = (x\, . . . , x r ) E C, and X = {ni, na, . . . , n a } C {1, 2, . . . , r} 
of size \X\ — a be given. Let the elements of X be enumerated in increasing order, 
i.e. 

1 $5 n\ < «2 < ■ ■ • < n a ^ r. 

Denote by u>x the vector 



Then we denote 



Wx — {Xni : ) • • * : x n a ) £ C • 



A x H d = A"" 1 (w x ) . 



We immediately obtain the following result. 



lr The coefficients {a^J may be readily obtained by the Cauchy residue formula 

a jt = 1 lim (±Y^- e \(*-^ 

3 ' (dj-l-ey.z^TjKdzJ { q„, w (z) 
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Lemma 6.6. For all w E C r and X C {1,2, ... ,r}, we have Ax (w) E V w . 
Moreover, letting a — \X\ we have 

SM(A x (w)) = A a - 1 (w x ) — = 1 - T ^. (6.9) 

Finally, (w,Ax (w)) is a continuous section ofV r . 

6.3. Constructing a basis. The following result is well-known, see e.g. [9, Propo- 
sition 35]. 

Theorem 6.7. Denote Nj = {1, 2, . . . ,j} for j = 1,2, ...,r. Then for every 
w E C r , the collection 

is a basis for V w . 

There are various proofs of this statement. Below we show how to construct sets 
which do not necessarily remain basis for all w E C r , but only for w in a small 
neighborhood of a given ido £ C r . Theorem 16. 71 will then follow as a special case of 
this construction. 

Informally, if two coordinates x\ and Xj can collide, then it is necessary to allow 
them to be glued by some element of the basis, i.e. we will need Ax (w) where 
i,j E X (in Theorem 16.71 all coordinates might be eventually glued into a single 
point because w is unrestricted.) In order to make this statement formal, let us 
introduce a notion of configuration, which is essentially a partition of the set of 
indices. 

Definition 6.8. A configuration C is a partition of the set N r — {1, 2, . . . ,r} into 
s = s (C) disjoint nonempty subsets 

U s i=1 Xi = N r , \X i \=d i >0. 

The multiplicity vector of C is 

T(C) = (d 1 ,...,d s ). 

Every configuration defines a continuous family of divided differences as follows. 

Definition 6.9. Let a configuration C = {Xj} s ^\ ■ Enumerate each Xj in increas- 
ing order of its elements 

Xj = [n{ < n\ < ...<.} 
and denote for every m = 1,2, ... ,dj 

Xj, m = f {n{ : k = 1, 2, . . . , m} . 
For every w E C r , the collection Be (w) C V w is defined as follows: 

Now we formally define when a partition is "good" with respect to a point w E C r . 

Definition 6.10. The point w = (x\, . . . ,x r ) E C r is subordinated to the configu- 
ration C — {Xj}j^l if whenever = xi for a pair of indices k 7^ I , then necessarily 
k, I E Xj for some Xj . 
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Now we are ready to formulate the main result of this section. 

Theorem 6.11. For a given wo € C and a configuration C, the collection Be (wq) 
is a basis for V Wo if and only if Wo is subordinated to C. In this case, Be (w) is a 
continuous family of bases for V w in a sufficiently small neighborhood of Wq ■ 

Let us first make a technical computation. 

Lemma 6.12. For a configuration C and a point w € C r , consider for every fixed 
j = 1, . . . , s (C) the set 

Sj = f {A Xjim (w)}^ =1 . (6.10) 

(1) Define for any pair of indices 1 ^5 k ^ £ ^ dj the index set 

Xj,k-.e = {"I < »;} Xj = X^-.di = x j4j . 

Then 

A x 3 , k .. e (w) e spaaSj. 

(2) For an arbitrary subset Y C Xj (and not necessarily containing segments 
of consecutive indices), we also have 

Ay (w) £ span Sj . 

Proof. For clarity, we denote yi — x n i and [k : £] = Axj k . e {w). By (|6.3p we have 
in all cases (including repeated nodes) 

(ye - Vk) [k : £} = [k + 1 : £} - [k : i - 1] . (6.11) 

The proof of the first statement is by backward induction on n = £ — k. We start 
from n = dj, and obviously [1 : dj] € Sj. In addition, by definition of Sj we have 
[1 : m] G Sj for all m = l,...,dj. Therefore, in order to obtain all [k : £} with 
I — k = n— 1, we apply (|6.11| several times as follows. 



[2 : n] = [y n - Vl ) [1 : n] + [1 : n - 1] 
[3:n+l] = {y n+l ~y 2 ) \2:n + l) + [2:n] 

[dj -n + 2: dj] = [y^ ~ y^-n+l) [dj -n+l : dj] + [dj - n + 1 : dj - 1] 

Here the symbol under a term means that the term is taken directly from the 
previous line, while — indicates that the induction hypothesis is used. In the end, 
the left-hand side terms are shown to belong to span Sj . 

In order to prove the second statement, we employ the first statement, (|6.8j) and 
Proposition EH Item[U □ 



Proof of Theorem \6.11\ In one direction, assume that wq = (x\, . . . ,x r ) is subor- 
dinated to C. It is sufficient to show that every element of the standard basis (|2.2p 
belongs to spanj^c (wq)}. 

Let Tj £ T(wq), let dj be the corresponding multiplicity, and let Yj C N r denote 
the index set of size dj 

Yj = {i : xi = Tj} . 
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By the definition of subordination, there exists an element in the partition of C, say 
Xk, for which Yj C Xk- By Lemma \6. 121 we conclude that for all subsets Z CYj, 

I X I 

A z (w ) € span { Ax feim (w )} m ^ C span{i3 c (w )} ■ 
By (|6.7p , (wo) is nothing else but 

A Z (w ) = A^ 1 f j = pjl-tfl'M) (, - TJ ) . 

\ x|Z| / 

This completes the proof of the necessity. In the other direction, assume by contra- 
diction that Xk — xi — t but nevertheless there exist two distinct elements of the 
partition C, say X a and Xp such that k <E X a and £ 6 Xp. Let the sets {Sj}^l be 
defined by (16. 10[) . Again, by Lemma \6 . 1 21 and (|6 . 7[) we conclude that 

5 (x — t) <E span S a n span S 1 ^ . 

But notice that Be (wq) = \Jj=l Sj and Y^j=i — therefore by counting 
dimensions we conclude that 

dim span {Be (wq)} < d, 
in contradiction to the assumption that Be (wo) is a basis. 

Finally, one can evidently choose a sufficiently small neighborhood U C C r of wq 
such that for all w € U, no new collisions are introduced, i.e. w is still subordinated 
to C. The continuity argument f Lemma I6.6[) finishes the proof. □ 

Remark 6.13. Another possible method of proof is to consider the algebra of ele- 
mentary fractions in the Stieltjes space S ri and use the correspondence (|6.9[) . 

As we mentioned, Theorem 16.71 follows as a corollary of Theorem 16.111 for the 
configuration C consisting of a single partition set N r . 

6.4. Resolution of collision singularities. Let /ij e S* C C 2r be given, and let 
(wo,go) € V r be a, solution to the (rank- restricted) Prony problem. The point wo is 
uniquely defined up to a permutation of the coordinates, so we just fix a particular 
permutation. Let T (wo) = (ri, . . . , t s ). 

Our goal is to solve the rank-restricted Prony problem for every input /1* € C 2r in a 
small neighborhood of /Iq. According to Theorem 15. 71 this amounts to a continuous 
representation of the solution (z) = g' t ^ = TM.^ (/•**) to the corresponding 
diagonal Padc approximation problem as an element of the bundle V T - 

Define S = minj=y |rj — Tj\ to be the "separation distance" between the clusters. 
Since the roots of Q M » depend continuously on fi* and the degree of Q^, does not 
drop, we can choose some /ij sufficiently close to jUq, for which 

(1) all the roots of Q M * (z) are distinct, and 

(2) these roots can be grouped into s clusters, such that each of the elements 
of the j-th cluster is at most (5/3 away from Tj. 



GEOMETRY AND SINGULARITIES OF THE PRONY MAPPING 



20 



Enumerate the roots of Q^* within each cluster in an arbitrary manner. This 
choice enables us to define locally (in a neighborhood of r algebraic functions 
X\ (/x*) , . . . , x r (/**), satisfying 

(z)=n(z- Xj o«*)). 

Then we extend these functions by analytic continuation according to the above 
formula into the entire neighborhood of /ig. Consequently, 

w(ff) = f (xi (/i*) ,...,x r (/i*)) 
is a continuous (multivalued) algebraic function in a neighborhood of /ig , satisfying 

w (hq) = w . 

After this "pre-processing" step, we can solve the rank-restricted Prony problem in 
this neighborhood of /ig, as follows. 



Algorithm 1 Solving rank-restricted Prony problem with collisions. 

Let /ig G S* C C 2r be given, and let (wo,5o) G V r be a solution to the 
(rank-restricted) Prony problem. Let wq be subordinated to some configuration C. 

The input to the problem is a measurement vector fx* = (mo, . . . , ?7i2r-i) G C 2r , 
which is in a small neighborhood of /ig. 

(1) Construct the function w — w (/i*) as described above. 

(2) Build the basis B c (w) = {A X]ll M}^^ for V w . 

(3) Find the coefficients {Pj,i}jZ^' "' s \c) sucn t\v&t 

by solving the linear system 
Pj,t M ^x.j. (w) (x k ) =m k (= f x k g (w) {x)dxj, k = 0, 1 



} * * • ) 



2r-l. 



(6.12) 



Theorem 6.14. XTie coordinates {/3j,i} of the solution to the rank-restricted Prony 
problem, given by Algorithm \6Ji\ are (multivalued) algebraic junctions, continuous 
in a neighborhood of the point /ig . 

Proof. Since the divided differences Ajj (w) are continuous in w, then clearly for 
each k = 0, 1, . . . , 2r — 1 the functions 

u jAk (w) = A jti (w) (x k ) = A 1 ' 1 (w Xj , e ) (x k ) 

are continuou^ in w, and hence continuous, as multivalued functions, in a neighbor- 
hood of fj,Q. Since Be (w (/i*)) remains a basis in a (possibly smaller) neighborhood 

2 In fact, Vjj^h (w) are symmetric polynomials in some of the coordinates of w. 
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of /Jq, the system (|6.12ll . taking the form 

^ Vj,t,k (w) /3j,t (w) = rrik, k = 0, 1, . . . , 2r - 1 

3,t 



remains non-degenerate in this neighborhood. We conclude that the coefficients 
{Pj,£ (w (/j,*))} are multivalued algebraic functions, continuous in a neighborhood 



In this section we shall restrict ourselves to the real case. Notice that in many 
applications only real Prony systems are used. On the other hand, considering the 
Prony problem over the real numbers significantly simplifies some constructions. In 
particular, we can easily avoid topological problems, related with the choice of the 
ordering of the points x\, . . . , Xd € C. So in a definition of the real Prony space RVd 
we assume that the coordinates x%, . . . ,Xd are taken with their natural ordering 
x\ < X2 < • • • < Xd- Accordingly, the real Prony space RVd is defined as the 
bundle (w,g), w G Y[d c e RV W . Here Y[d * s the prism in M. d defined by the 

inequalities X\ < X2 < ■ ■ • < Xd, and RV W is the space of linear combinations with 
real coefficients of ^-functions and their derivatives with the support {x\, . . . ,Xd}, 
as in Definition 12.41 The Prony, Stieltjes and Taylor maps are the restrictions to 
the real case of the complex maps defined above. 

In this paper we just point out a remarkable connection of the real Prony space 
and mapping with hyperbolic polynomials, and Vieta and Vandermondc mappings 
studied in Singularity Theory (see [TJ [131 HH [H] an d references therein). 

Hyperbolic polynomials (in one variable) are real polynomials Q(z) — z d +'Y^ d = i ^jZ d 
with all d of their roots real. We denote by Td the space of the coefficients 
A = (Ai,...,Ad) C M d of all the hyperbolic polynomials, and by Td the set of 
A G Td with Ai = 0, IA2I < 1. Recalling (|2.3p . it is evident that all hyperbolic 
polynomials appear as the denominators of the irreducible fractions in the image 
of RVd by SM. This shows, in particular, that the geometry of the boundary dT 
of the hyperbolicity domain T is important in the study of the real Prony map 
VM.: it is mapped by VM. to the boundary of the solvability domain of the real 
Prony problem. This geometry has been studied in a number of publications, from 
the middle of 1980s. In [13] V. P. Kostov has shown that T possesses the Whit- 
ney property: there is a constant C such that any two points Ai,A2 G T can be 
connected by a curve inside T of the length at most C|| A2 — Ai||. "Vieta mapping" 
which associates to the nodes x\ < X2 < • • • < Xd the coefficients of Q(z) having 
these nodes as the roots, is also studied in [13) . In our notations, Vieta mapping is 
the composition of the Stieltjes mapping SA4 with the projection to the coefficients 
of the denominator. 

In [T] V.I. Arnold introduced and studied the notion of maximal hyperbolic poly- 
nomial, relevant in description of T. Furthermore, the Vandermonde mapping 
V : R d -> R d was defined there by 



of fj*. 



□ 



7. Real Prony space and hyperbolic polynomials 
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with at, . . . , ad fixed. In our notations V is the restriction of the Prony mapping 
to the pairs (w, g) € RVd with the coefficients of g in the standard basis of RV W 
fixed. It was shown in [I] that for at, ■ ■ • ,a-d > V is a one-to-one mapping of 
Y[d to its image. In other words, the first d moments uniquely define the nodes 
x\ < X2 < ■ ■ ■ < Xd- For at, • ■ ■ ,a>d with varying signs, this is no longer true in 
general. This result is applied in [1] to the study of the colliding configurations. 

Next, the "Vandermonde varieties" are studied in [T], which are defined by the 
equations 

atxt + . . . + a d Xd = a.1, 

aix[ + . . . + a d x p d = a e . 

It is shown that for a\, . . . ,a d > the intersections of such varieties with FJ d are 
either contractible or empty. Finally, the critical points of the next Vandermonde 
equation on the Vandermond variety are studied in detail, and on this base a new 
proof of Kostov's theorem is given. 

We believe that the results of [TJ [T3] and their continuation in |14l [T5] and other 
publications are important for the study of the Prony problem over the reals, and 
we plan to present some results in this direction separately. 



Appendix A. Proof of Theorem 13.51 

Recall that we are interested in finding conditions for which the Taylor mapping 
TA4 : S d —> T d is invertible. In other words, given 

2d-l (\\ k + x 
k=0 

we are looking for a rational function R (z) S S d such that 

Write R(z) = gf} with Q (z) = £^ =0 cy^ and P (z) = £)£o Multiplying 
(|A.1[) by Q (z), we obtain 

Q(*)S(*)-P(,) = -£ r + -g__ + .... (A.2) 

Proposition A.l. The identity (|A.2[1 . considered as an equation on P and Q with 
deg P < deg Q < d, always has a solution. 

Proof. Substituting the expressions for S, P and Q into (|A.2[) we get 

(co + ctz + ■ ■ ■ + c d z d ) + ^ + . . . ) -bo bd-^- 1 = . . . (A.3) 

The highest degree of z in the left hand side of (jA.3|) is d — 1 . So equating to zero 
the coefficients of z s in (|A.3[) for s = d — 1, . . . , — d we get the following systems of 
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equations: 





m 



m 
m 1 





Cl 




bd-1 




C2 




bd-2 








. bo _ 



_m mi ... m<j-i_ 
From this point on, the equations become homogeneous: 



(PEN 



mo 


mi 


m d 




co 




"0" 


mi 


m 2 . . 


m d+ i 




Cl 







m d -i 


m d . . 


■ m 2 d-i_ 




, C d_ 








The homogeneous system (|AI**|) has the Hankel-type d x (d + 1) matrix Md = 
(TOi+j) with ^ i ^ e? — 1 and ^ j ' ^ <i. This system has d equations and e? + 1 
unknowns Co, . . . , c<j. Consequently, it always has a nonzero solution cq, . . . , C4. Now 
substituting these coefficients cq, of Q into the equations (|Al^r|) we find the 

coefficients &o, • ■ • , bd~i of the polynomial P, satisfying (|Al*p , Notice that if Cj = 
for j ^ £ + 1 then it follows from the structure of the equations (|AI*j) that 6j = for 
j > £. Hence these P, Q provide a solution of (|A.2j) , satisfying deg P < deg Q < d, 
and hence belonging to Sd- □ 



However, in general (|A.2|) does not imply (|A.1|> . This implication holds only if 
deg Q = d. The following proposition describes a possible "lost of accuracy" as we 
return from (|A.2|1 to (|A.1|) and deg Q < d: 

Proposition A. 2. Let (|A.2[) be satisfied with the highest nonzero coefficient of Q 
being C£, £ < d. Then 

P{z) dt d 2 



S{z) 



Q{z) 



yd+l+2 



(A.4) 



Proof. We notice that if the leading nonzero coefficient of Q is eg then we have 

7; = — ( 1 Q_x j ) = — (/0 + fl~ + ■■■)■ 

Q z* Ci + + . . . z l z 
So multiplying (|A.2|) by ^ we get (|A.4[) . □ 



Proof of Theorem \3.5\ Assume that the rank of Md is r < d, and that \M r \ 7^ 0. 
Let us find a polynomial Q(z) of degree r of the form Q(z) = z r + X)j=o c i z "'' 
whose coefficients satisfy system (|AI**j) . Put c r = (co, . . . , c r -i, 1) T and consider a 
linear system M r c r = 0. Since by assumptions |Af r | ^ 0, this system has a unique 
solution. Extend this solution by zeroes, i.e. put = (cq, . . . , c r _i, 1,0,..., 0) T . 
We want to satisfy (|AI^| . which is MdCd = 0. This fact is immediate for the 
first r rows of Md- But since the rank of Md is r by the assumption, its other rows 
are linear combinations of the first r ones. Hence satisfies (|A1**[) . 

Now the equations (|A1*[) produce a polynomial P(z) of degree at most r — 1. So 
we get a rational function -R(z) = 1^4 € <S r C 5^ which solves the Pade problem 
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(TO) , with degQ(z) = r. Write R(z) = Y,T= a k(^) k+1 ■ % Proposition [O we 
have rrik = otk till k = d + r — 1. 

Now, the Taylor coefficients a k of -R(-z) satisfy a linear recurrence relation 



2~^c s m fe _ s , k = r,r + l, (A. 5) 



Considering the rows of the system M^A = we see that m k satisfy the same 
recurrence relation (|A.5I) till fc = d + r — 1 (we already know that m k — a k till 
k = d + r — 1). We shall show that in fact m k satisfy (|A.5[) till fc = 2d — 1. 

Consider a d x r matrix formed by the first r columns of Md, and denote its 
row vectors by Vj = (m,-,o, • ■ • , mj )r _i), i = 1, . . . , d — 1. The vectors Vj satisfy 

r 

Vj = -^c s Vj_ s , i = r, ...,d-l, (A. 6) 

s=l 

since their coordinates satisfy (|A.5[) till fc = d + r — 1. Now vo, . . . , v r _i are linearly 
independent, and hence each v,, i = r, . . . , d — 1, can be expressed as 

r-l 

Vj = ^7i, s v 5 .. (A. 7) 

s=0 

Denote by Vj = (m^o, ■ • ■ , i = 1, . . . , d — 1 the row vectors of M^. Since by 

assumptions the rank of Md is r, the vectors Vj can be expressed through the first 
r of them exactly in the same form as v^: 

r-l 

Vi = ^7 M v s , i = r, ...,d-l. (A. 8) 

s=0 

Now the property of a system of vectors to satisfy the linear recurrence relation 
(|A.6[) depends only on the coefficients 7$ g in their representation (jA.7[) or (|A.8I) . 
Hence from (|A.6|) we conclude that the full rows Vj of Md satisfy the same recurrence 
relation. Coordinate- wise this implies that m k satisfy (|A.5|) till k = 2d — 1, and 
hence m k = a k till fc = 2d — 1. So R(z) solves the original Problem 13. II 

In the opposite direction, assume that R(z) solves Problem 13. 1[ and that the rep- 
resentation R(z) = e S r C Sd is irreducible, i.e. degQ = r. Write Q(z) — 
z r +Y2j=o c j z ~ ' ■ Then m k , being the Taylor coefficients of R(z) till k — 2d— 1, satisfy 
a linear recurrence relation (IA.5|) : m k — ~ 53s=l c s TO fc-s, k = r, r + 1, . . . , 2d — 1. 
Applying this relation coordinate-wise to the rows of Md we conclude that all the 
rows can be linearly expressed through the first r ones. So the rank of Md is at 
most r. 

It remains to show that the left upper minor \M r \ is non-zero, and hence the rank 
of Md is exactly r. 

By Proposition 13.31 if the decomposition of R(z) in the standard basis is 

• * (-1)^(1-1)! 

R ( z ) = 2^ 2^ a ^-! — , _ re — - 

i=xi=x y z x j> 

where X^=i = r & n d are pairwise distinct, then the Taylor coefficients of 
R(z) are given by (|1.5p . Clearly, we must have aj^-i 7^ for all j = 1, . . . , s, 
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otherwise deg Q < r, a contradiction. Now consider the following well-known rep- 
resentation of M r as a product of three matrices (see e.g. [7]): 

M r = V(x 1 ,d 1 ,...,x s ,d s ) x diag{A,}* =1 x V (xi, d 1: . . . , x s , d s ) T , (A.9) 

where V (. . . ) is the confluent Vandermonde matrix (|4.1[) and each Aj is the fol- 
lowing dj x dj block: 



A 3 ^ 



a j,0 a J,l ' ' ' ' ' ' a j,dj- 






(VK^-i ••• 
.a^dj-i ••• ••• 

The formula (|A.9[) can be checked by direct computation. Since {xj} are pairwise 
distinct and dj^-i for all j = 1, . . . , s, we immediately conclude that \M r \ ^ 0. 
This finishes the proof of Theorem 13.51 □ 
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